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Abstract 

We parameterize the initial flux composition of high energy astrophysical neutrinos 
as (<^g : <I>J^ : = (1 : n : 0), where n characterizes the source. All usually assumed 
neutrino sources appear as limits of this simple parametrization. We investigate how 
precise neutrino telescopes can pin down the value of n. We furthermore show that 
there is a neutrino mixing scenario in which the ratio of muon neutrinos to the other 
neutrinos takes a constant value regardless of the initial flux composition. This occurs 
when the muon neutrino survival probability takes its minimal allowed value. The 
phenomenological consequences of this very predictive neutrino mixing scenario are 
given. 
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1 Introduction 



Neutrino mass and mixing influences many phenomenological aspects of particle and as- 
troparticle physics. In particular, any observed flavor mix of neutrinos Ve-, i^/i and Ur is 
usually not the original flavor mix, but rather a modified one. This is a consequence of 
the non-trivial structure of the leptonic mixing, or the Pontecorvo-Maki-Nakagawa-Saki 
(PMNS) mixing matrix. In this letter we wish to focus on the interplay of neutrino mix- 
ing and the flavor composition of ultra high energy (UHE) astrophysical neutrinos, whose 
detection is the goal of neutrino telescopes such as IceCube f], currently under construc- 
tion, or the planned KM3Net facility [2]. Many recent works discussed aspects of this 
problem [3142T]. an overview can be found in [22] . 

The composition of the flux detectable in neutrino telescopes depends on two things - 
(i) the initial flavor composition at the source, and (ii) the standard (and possibly non- 
standard) parameters of neutrino mixing. This means that there are two types of studies 
one may pursue with neutrino telescopes: 

(a) Neutrino physics: the extraction of neutrino properties from measurements of flavor 



(b) Astrophysics: the identification of the initial flavor composition of the flux of UHE 
neutrinos and using that to probe the type of source and its properties. 

A large plethora of literature exists for item (a), where authors have studied the potential of 
probing the neutrino mixing angles using the flavor ratios. In this paper we will concentrate 
on item (b) and expound the potential of UHE neutrino measurements to decipher the 
flavor composition of the UHE neutrinos at their source (see Ref. [20] for related recent 
analyses). This could lead to a better understanding of the properties of the astrophysical 
source responsible for the production of these neutrinos. As will be discussed in detail in 
Section 3, different astrophysical situations could give rise to a variety of flavor ratios of 
the neutrinos. Since the usually assumed high energy neutrino sources do not generate tau 
neutrinos in any appreciable amounts, we neglect their presence in the initial UHE neutrino 
flux and propose a simple one parameter parameterization of the initial flux composition 
of the neutrinos: 



All discussed sources such as pion, muon-damped, charm or neutron beams are simple limits 
of this parameterization. With such a parameterization the experimental determination of 
the initial flux composition is simply an extraction of the parameter n from the observed 
flux composition on Earth. We will define two kinds of measurable flux ratios and calculate 
the predicted values of these ratios at Earth for all values of n from — oo, a range 
which covers all possible UHE neutrino flux sources. We show that the uncertainty on the 
predicted values of these flux ratios due to their dependence on the oscillation parameters 
threatens to wash out any sensitivity of the neutrino telescopes to n. We define a very 
simple function and show quantitatively the ranges of n which could be ruled out by 
the data from UHE observations. 
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We also point out, for the first time, a special mixing scenario for which the observable 
ratio of muon neutrinos to the other neutrinos is independent of the initial mixing scenario. 
This occurs when the averaged survival probability for muon neutrinos takes its minimal 
allowed value of |. It illustrates nicely that the ratio of muon neutrinos with the other 
fiavors alone may not be a good discriminator of different sources. 

Our paper is build up as follows: in Section |2] we summarize the neutrino mixing phe- 
nomenology framework at neutrino telescopes. Section [3] sees the discussion of several 
neutrino sources whose common features lead to the parameterization in Eq. ([T]). The 
peculiar mixing scenario with a constant muon neutrino fiux with respect to the total fiux 
is presented in Section |H The experimental determination of n is investigated in Section 
[5l before we sum up and conclude in Section [61 



2 Neutrino Mixing and Neutrino Telescopes 

Astrophysical sources will generate fiuxes of electron, muon and tau neutrinos, denoted by 
$g, $° and respectively. As a consequence of non-trivial neutrino mixing, it is not this 
initial fiux composition which arrives at terrestrial detectors. In fact, what is measurable 
is given b}0 



(2) 



( $e ^ 




P 

^ e/x 


Per \ 


1 






P 

MM 










P 


Prr 1 





where the neutrino mixing probability is 

PaP = Pfia 



(3) 



and U is the lepton mixing matrix. The current best-fit values as well as the allowed la, 
2cr and So" ranges of the oscillation parameters are [23] 
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0.312 0.016 0.466 

0.294^0.331 0.006 ^0.026 0.408^0.539 

0.278^0.352 0.000 ^0.036 0.366^0.602 

0.263^0.375 0.000 ^0.046 0.331 ^0.644 



These mixing angles can be related to elements of the PMNS mixing matrix via 
/ C12C13 S12C13 sise""^ \ 



U 



-S12 C23 - C12 S23 Si3 e^^ C12 C23 - S12 S23 Si3 e'^ S23 Ci3 
\ ■S12 S23 — C12 C23 Si3 e**^ — C12 S23 — S12 C23 Si3 e*'' C23 Ci3 y 



(4) 



(5) 



^Everywhere in this paper we denote the fiuxes at the source with a superscript ($'^), while the 
observed fluxes will be denoted without any superscript. 
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where Cij 



cos 6ij, Sij = sin 6ij. The CP phase S is unknown. Because the travelled 



distance of the neutrinos is much larger than the oscillation length An E/Am'^, the mass- 
squared differences Am? drop out of the mixing probabilities. Furthermore, solar neutrino 
mixing is neither maximal, nor zero or 7r/2, and due to these reasons no transition prob- 
ability Pai3 with a P is zero and no survival probability Paa is one [iTj. Consequently, 
high energy astrophysical neutrinos will always mix. To be precise, at la and 3a the entries 
of the flavor conversion matrix P^,^ are 

/ 0.529 ^0.578 0.178^0.296 0.158 ^0.275 
0.341^0.443 0.354^0.394 
0.354^0.469 

/ 0.486^0.612 0.127^0.344 0.118 ^0.335 
0.333^0.508 0.304^0.403 
0.333 0.525 
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(at la) , 
(at 3a) . 



(6) 



0.333 and P^r = 0.333 



The consequences and phenomenology of the special values P^^ 
will be dealt with in Section HI 

As is obvious from the ranges of the mixing parameters given above, a good zeroth order 
description of neutrino mixing is 
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Therefore, it proves in particular useful to expand in terms of 



f/psl and e 
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e + 0{e^). The result of the expansion for the flavor mixing matrix 
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where the universal first [HIITT] and second [l7l[T8] order correction terms are 



A = - sin^26'i2e + 7 sin 46112 cos6\Ue3\ — - (V2 cos5 If/esl 



4e 



(cos 26^12 e — sin 26*12 cos 5|f/e3|)^ ~3e2 + (2v^ cos5|f/e3| -e) 



(10) 



We have also given the expressions for the value sin^ ^12 



|. Note that A is positive 
semi-definite. The same holds for the third expansion parameter 
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which is usually only a sub-leading correction. The parameter A can seldomly be larger 
than the first order term A. To be quantitative, 

at la : -0.046 < A < 0.069 , A^ < 0.060 , 0.003 < f/^ < 0.015 , 

" " -2~ ~~ (12) 

at 3a: -0.101 < A < 0.112 , A < 0.162 , [/^ < 0.028 . 

The dependence of the expansion parameters on Qyi is weak. From the expressions for A 

2 

and A it is clear that their dependence on the atmospheric neutrino mixing angle 6*23 is 
stronger than the one on \Ue-i\ cos 5 [T9] . 



3 Initial Flux Compositions and their Parametriza- 
tion 

We now turn to the observable flavor ratios [2ll[25]. The most frequently considered and 
experimentally most accessible is the ratio of flux of muon neutrinos to that of all other 
flavors: 

^^ = -^- (13) 

Often one considers also the ratio ($e + which is simply T j (1 — T). Whereas muon 
neutrinos with their characteristic tracks (continuous loss of energy via Cerenkov radiation) 
can be somewhat easily distinguished from electron and tau neutrinos, the distinction of 
the latter two is more difficult. If possible (e.g., via electromagnetic vs. hadronic showers, 
or tau-induced lollipop/double bang events), one may consider their ratio as well: 

(14) 

We will use these ratios T and R in this paper. When neutrinos and anti-neutrinos are 
not distinguished, then two independent ratios are sufficient to fully determine the flavor 
content. For instance, if both T and R are known, the ratio of muon neutrinos to tau 
neutrinos 5* is related to T and R via T = 5'/(S' + i? + 1). 

In what regards initial flux compositions, typically one considers pionic beam-dump-like 
sources, which have an initial flux composition of 

($0 : $J : $0) ^ (1 ^ 2 : 0) . (15) 

They arise because hadronic processes in a cosmic ray source produce pions. If the medium 
in such a source is opaque to muons, then one speaks of muon-damped cases [26], which 
have a composition of: 

($" : $2 : $?) = (0 : 1 : 0) . (16) 

Another possibility are so-called neutron beams, which result when photo-disintegration 
creates neutrons from nuclei [27], and therefore 

: $J : $°) = (1 : : 0) . (17) 
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Finally, at high energies semileptonic decays of charm quark^ may generate an initial 
flavor ratio of [28] 

(<|.° : $0 : $0) = (1 : 1 : 0) . (18) 

We see that in all sources no tau neutrinos are generated initially. There is in fact always a 
small component ("prompt z/,-") from decays such as Dg ^ t u^-, however this contribution 
is at most at the permille level and can safely be neglected. We stress however that all of 
the above sources can be expected to be "impure", i.e., there will be small deviations from 
the idealized compositions given above [9 t ll3 p i8]. As one example, Ref. [13] has argued that 
because of the wrong helicity polarization of the muons in pion decay the z/^ have actually 
a softer spectrum, and thus the effective count of muon neutrinos is reduced. This depends 
on the injection spectrum, which is described by the power law For the canonical 



value a = 2 it was found that effectively ($^ : $Ji : $^) = (1 : 1.86 : 0). It was shown [18] 
that inclusion of leptonic and semi-leptonic kaon decays, as well as of heavy flavor decays 
does not change this ratio appreciably. 

One may also consider GZK- neutrinos [29j, whose flavor content changes with energy. Be- 
low about 100 PeV (10*^ GeV) one has (1:0:0), whereas for higher energies (1:2:0). 

Given all the above discussion, we are lead to parameterize the initial flux composition of 
high energy neutrinos simply as 

($0 ..^o^:^0) = (l:n:0). (19) 

The experimental determination of the initial flux composition is therefore simply an ex- 
traction of the parameter n from the observed flux composition on Eartlil. From Eq. ([19]) 
one can obtain (pure) pion sources when n = 2, neutron beams for n = 0, charm sources 
for n = 1 and muon-damped sources when one takes the limit n — )■ cxd. 
Let us attempt to express the flux ratios T and i? in a source-independent way. We take 
advantage of the parametrization in Eq. ( fT9l) and use the expansion of the probabilities 
Pap in Eq. ([9]). First we consider the ratio of muon to all neutrinos and flnd 

$ P +nP 



T 



$tot ~ 1 + ^ 



1 + n 



■ ' 1 - - I c?2 + A (1 - n) + - A + f/M 1 - - 



(20) 



.2 V 2/ ' '2 V 2, 

We note the following points from this expression: 

• for n = 2, which corresponds to pionic sources, there is no explicit dependence on 
9i2 and also corrections due to U"^ vanish. The uncertainty on T is hence expected 
to be amongst the lowest for this case; 

• for n = 1, which corresponds to charm sources, there is no flrst order correction A, 
and thus corrections to the zeroth order expression are only of quadratic order in the 
small expansion parameters 6*23 — 7r/4 and |f/e3|; 



^For bottom quark decays the production rate is suppressed by one order of magnitude. 
3Ref. [6] has proposed ($° : $° : $°) = (sin^ cos^ C : cos^ ^ cos^ C : sin^ C)- 
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la 


3cr 


la 


3(T 


pion 
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0.324^0.355 


0.323^0.387 


0.893 ^ 1.247 


0.820 ^ 1.450 


charm 


1 


0.303 -^0.323 


0.298 ^ 0.348 


1.121 ^ 1.581 


0.998 ^1.853 


neutron 





0.178^0.296 


0.127^0.345 


1.922 ^3.618 


1.449 ^5.090 


muon-damped 


oo 


0.341 ^ 0.443 


0.333^0.508 


0.462 ^0.814 


0.344^1.071 



Table 1: Ranges of the flux ratios T = $^/$tot and R = $e/'^'r for an initial flux compo- 
sition of ($g : : = (1 : n : 0). The case n = 2 correspond to pion sources, n = 1 to 
charm sources, n = to neutron beams and n — )■ oo to muon-damped sources. We have 
inserted the la and 3a ranges of the oscillation parameters. 



In the zeroth order approximation, which corresponds to /i-r symmetry yielding Oi^ = 
and 623 = such that the correction terms A, A and If are all vanishing, Eq. f EU]) reduces 
to 

2 (pion source) 
1 (charm source) 
(neutron beam) 
00 (muon-damped) 

where we have also given the simple expressions for the special values of n corresponding 
to the usual sources one considers. Hence, we note that all sources have similar predicted 
T around 0.3. Indeed, in Section |4] we will present an extreme scenario in which T = | 
independent of n. We further note that the predicted T is highest for n = 00, lowest for 
n = and intermediate for n = I and 2. In general, T slightly increases with n. 
Similarly, the ratio R of electron and tau neutrinos is 

~ Per + nP^r 

l-2c22S?2(l-f)+nA-2^/2 (1-1) (21) 
^t + cL^L (l-f)-A-f AVf/2 (l-f)- 

Here the value n = 2 removes again several correction terms, in particular the terms with 
explicit dependence on 612, which leads to small uncertainty. For the case of /i-r symmetry 
one flnds 

2 (pion source) 
1 (charm source) 
(neutron beam) 
00 (muon-damped) 
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The spread for the different values of n is larger than for T. For instance, if n = 2 then 
i? = 1 at leading order and between 0.820 -v- 1.450 for the 3a range, while n = 1 leads to 
= YY at leading order and between 0.998 ^ 1.853 for the 3cr range. In general, R slightly 
decreases with n. 

In Fig. [T]we show T and i? as a function of n. In Table [1] we present the la and 3a ranges 
of predicted T and R for the 4 benchmark cases of n corresponding to pion, charm, neutron 
and muon-damped sources. The la and 3a ranges of T and R correspond to current la 
and 3a ranges of the oscillation parameters. In general one notes again that T is not an 
ideal discriminator in order to distinguish the neutrino source, because it is rather close 
to I for all n. Therefore, at the present stage we can already suspect that the task of 
identifying n will be easier when R is added to the statistics. A more detailed analysis will 
be presented in Section 5. 



4 An Extreme Case: Minimal Survival Probabilities 

One observes from the allowed ranges of the oscillation probabilities in Eq. ([6]) that P^^ and 
Prr can take the value 0.333 when the 3a ranges of the oscillation parameters are inserted. 
It is worth noting that | is the minimal value that an averaged survival probability Paa 
can take if there are 3 fermion families. In general, in the presence of k flavors we have 

= I ■ (22) 

The minimum is obtained if and only if If/^jp = 1/k for all i, i.e., when all mixing matrix 
elements of a row take on the same valu^. 

Considering the case of A; = 3, let us discuss first the phenomenological consequences of 



Wfii\ = Y 1/3- The case of \Uri\ = yl/3 will be discussed at the end of this Section. First 
we note that unitarity of the PMNS matrix links the three entries of a row, which implies 



that only two independent constraints can result. The first one stems from \U^s\ = \/l/3 
and is 

sin^ e,s = I, I, ^ ^ (l + \Ue3\') . (23) 



31 - f/, 



e3| 



Inserting this in = Jl/3 gives the second independent constraint 



|2 



cos5tan2^i2 = ^ 2 [^esi ^ ^ - / ^_ ^ r [[/^g] 1 (24) 



f/e3K/2-3|t/e3| 



|2 




We immediately see that \Ue3\ should be rather large, sin^ 6*23 significantly less than | and 
cos 6 should be different from zero, or sin 6 different from 1. We display the phenomenology 
resulting from Eqs. ( 123| [2^ in Fig. |2l While sin^ 6'23 lies close to its lower 3a bound, sin^ 612 



Note that this concerns the elements of a row. Commonly it is assumed the elements of a column are 
identical (as for tri-bimaximal, or more general tri- maximal mixing [30 ). 
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Figure 1: Range of allowed values of T = $^/$tot (upper plot) and R = ^e/^r (lower 
plot) as a function of the flux composition parameter n, when the oscillation parameters 
are varied within their 3a ranges. 
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can approach its upper Icr bound from above, in which case |?7e3| is close to its current 
upper bound and 6 less than one (or larger than 5). Vanishing ^13 leads to maximal 612 
and is disallowed. 

The condition If/^^jp = | together with unitarity of the PMNS matrix fix the form of the 
matrix It/mp, namely 



\Ue3 

\Uai\ = I 3 3 3 I ■ (2^) 

Recall that the |t/Q,jp, and therefore the averaged probabilities Pap, depend on four inde- 




2 



pendent parameters. The condition |?7^j| = yl/3 fixes two of the four parameters and we 
have chosen |f/eiP and |f/e3p as the remaining two free and independent parameters. It is 
easy to see that from Eq. ( l25l) the relation 

Pfia = -Pa/. = ^ for a = e,fi,T (26) 

follows. With this relation one obtains a remarkably simple result for the ratio T as a 
function of n, namely 

rp _ _ -Pe/t + n _ 1 -^^N 

$tot 1 + n 3' ^ ' 

regardless of n! The origin of this simple result is of course that all neutrinos Vi have 
an equal share in i/^, i.e., their muon content is equally distributeclf]. Concerning i?, the 
general result as a function of n is rather lengthy. The allowed ranges of R are however 
significantly less than for the unconstrained case. For instance, R for pion sources lies 
between 1.23 and 1.40, to be compared with the general allowed range between 0.82 and 
1.45. An interesting case occurs for muon-damped sources, for which we find that R= 1. 

If the third row of U has identical entries, and therefore P^t = |, then 



Sin^ ^23 = ^ (2 -me,?)^l{l-\ |f/e3p ) (2^ 



and 

1 — 2 If/ P 

cos5tan20i2 = ' - ^ I t7t-\ " 7 1 • (29) 



t/e3U/2-3|f/e3| 
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While the second constraint is very similar to the one for \U^i\^ = |, the first condition 
implies here that sin^ 6^23 lies slightly outside its allowed 3a range. The value P™™ = 0.333 
found in Eq. is therefore only numerically close to |. In any case, since the identification 
of muon neutrinos is much simpler than of tau neutrinos, it makes little sense to discuss 
the scenario with a saturated Pi: 



jmm 

TT ' 



^Similarly, the ratio of muon neutrinos to the sum of electron and tau neutrinos is $^/($e + = i. 
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Figure 2: Phenomenology of a PMNS matrix with all elements of the second row equal 
to each other. Shown are the atmospheric neutrino parameter sin^^23 against If/gsl, the 
solar neutrino parameter sin^ 9i2 against 6 fi J different values of | Ue3 \ and sin^ 9i2 against 
sin^6'23. Also given are the current best-fit value and the la as well as 3a ranges of the 
oscillation parameters. 
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Figure 3: cis a function of Ufu for ntme = 2, for different values of errors on the flux 
ratios T (left panel) and R (right panel). All oscillation parameters are fixed in the fit. 



5 Experimental Distinction of Sources 

In this Section we explore the potential of the neutrino telescopes to ascertain the initial 
flux composition at the source. We stress that all results to be shown have been obtained 
by an exact calculation of the oscillation probabilities. We begin by returning to Fig. [T|, in 
which we display the current 3cr predicted ranges of T and R respectively, as a function 
of the flux composition parameter n. For each value of n we calculated the full range of 
predicted values of T and R by allowing the oscillation parameters to vary in their current 
3a allowed range [23]. The range of values of T and R for the different limiting cases are 
given in Table [TJ As discussed before, we note from the flgure that the T ranges for all 
n are overlapping. We also note that the spread in T is minimized for values of n close 
to 2, as for this case the explicit dependence on 612 vanishes (cf. Eq. fl20|) ). For the same 
reason, the spread in R is also minimized for values of n close to 2 (cf. Eq. fl2Tl) ). The 
spread in T for small n is somewhat larger than for large n. The spread in R for small n 
is seen to be signiflcantly larger than for large n. The large spread for small n of both T 
and R implies that neutron sources have the largest dependence on the neutrino mixing 
parameters. Note that while the possible spread induced by neutrino mixing is smallest 
for n between 1 and 2, it does not necessarily imply that the sensitivity to initial flux 
composition will be best for these cases. We also note that R appears to be a much better 
discriminator for n compared to T and hence inclusion of R in addition to T will make 
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the identification of tlie source easier. Tlie individual lines embedded in the band in Fig. [T] 
show the predicted T and R for several benchmark mixing scenarios, in which the values 
of the oscillation parameters are 

"TBM" the tri-bimaximal scenario: the oscillation parameters are s\v? 612 = |, sin^ ^23 = |, 
\Ue3\^ = and 5 = 0; 



: the oscillation parameters are sin 612 = 0.27, sin^ $23 = 0.35, 
= 0; 

: the oscillation parameters are sin 612 = 0.27, sin^ 023 = 0.65, 
= 0; 

"SPL" special scenario: for the special scenario discussed in the previous Section, we take 
sin^^ia = 0.352, sin^ ^23 = 0.348, jf/esp = 0.043 and 6 = 0. 

One can immediately see from Fig. [T]that as predicted in Section HI the special case yields 
T = 1/3 for all values of n. Therefore, if this case was the true mixing scenario chosen by 
Nature, then it would be impossible to conclude anything about n from T measurement 
at neutrino telescopes. One can also note that T for the EXl mixing case is very close to 
the SPL case and hence it would also be difficult to say anything about the source, if this 
set of mixing parameters turn out to be the true mixing angles. 
We next define a very simple function a^ 



"EXl" extreme scenario 1 
|f/e3p = 0.04 and 5 

"EX2" extreme scenario 2 
= 0.0 and 5 = 



2 ( Foata FTheory 



(30) 



where F can be either T or R and ^^Data" and ^^Theory" refer to observed and predicted 
fiux ratios respectively, and ap is the la experimental error on the relevant fiux ratio. We 
show results for fixed experimental errors of 20%, 10% and 5% for T and 30%, 20% and 
10% for R. We generate the '"F^iata" at certain "true" values n, denoted as Utrue, and at 
one of the benchmark mixing scenarios given above. This Foata is then fitted with FTheory 
corresponding to Ufu. The oscillation parameters in the fit are generally allowed to vary 
freely in the fit, apart from Fig. [3l 

In Fig. [3] we show the obtained as a function of nju for ntme = 2 and TBM mixing. 
All oscillation parameters are kept fixed at their TBM values in the fit for this figure. We 
reiterate that rifrue is defined as the value of n for which the Tdq^ (RData) are generated 
and Ufit is the value of n in Txheory {Rrheory)- The best-fit comes at njit = 2, as expected. 
We can see that extremely good sensitivity to n comes for AT/T = 5%, while reasonable 

^We are aware that this definition of is not strictly valid for a ratio as its error will not be Gaussian. 
Nonetheless we use it here for the sake of simplicity as the purpose of this paper is to roughly illustrate the 
potential of the neutrino telescope rather than to show exact numerical results for which one should work 
with the actual number of events calculated using a detailed code for the detector including threshold, 
efficiencies and so on. 
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M oscillatioD parameters free 




Figure 4: as a function of nju for ntme = 2, for different values of errors on the flux 
ratios R. All oscillation parameters allowed to vary freely within their current Sex limits 
in the fit. Corresponding for T is found to be zero for all values of Ufu and all errors 
AT/T. 



sensitivity to n is expected if AT/T = 10%. For AT/T = 20% one finds hardly any sensi- 
tivity at all, even in this ideal case, where the effect of oscillation parameter uncertainties 
have been neglected. If R could be measured at the neutrino telescope, then we expect 
good sensitivity even if AR/R = 30%. 

Now we allow the oscillation parameters to vary freely in the fit - the only restriction being 
that they are not allowed to take values beyond their current 3a limits. We find that, as 
expected from Fig. [U there is absolutely no n sensitivity if T is used. The reason, as 
discussed before, is the following: for TBM mixing parameters T^ata = 1/3 for ritrue = 2. 
This value of T can be reproduced by all values of n, as long as we are allowed to pick the 
most suitable set of oscillation parameters which lie within the current 3a limit. Therefore, 
for this case the TTheory always exactly reproduces the Toata and thus = for all nfn. 
For R the situation is slightly better. Roata = 1 for TBM mixing and ntme = 2. And 
hence, as predicted from Fig. [H values of nfn < 1 can still be disfavored. The reason being 
that for Ufit < 1, it is impossible to get the Rxheory close to 1. However, for nfn > 1, one 
can always find a set of oscillation parameters within the still allowed region which gives 
the same Roata- Fig. H] shows the result of the fit. 
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Marginalized over oscillation parameters| 





Marginalized over oscillation parameters 





Figure 5: For each ntme on the x-axis, the shaded region shows the values of nf^ which 
can be excluded at 2a. We assumed 10% (black/darkest), 20% (green/dark) and 30% 
(cyan/lightest) uncertainty on AR/R. The four panels are for the four benchmark values 
of the oscillation parameters assumed in the data. 



So far we have shown all sensitivity results assuming that we had ritrue = 2. We next 
show 2cr contour plots in the ntme — ^fit plane in Figs. [5] and [61 The way these figures are 
to be interpreted is the following: for each ntme on the x-axis, the shaded region shows 
the values of nja which can be excluded at 2a. Fig. [5] shows the sensitivity on n of the 
neutrino telescope using the more powerful flavor ratio R, while Fig. |6] gives the corre- 
sponding reach by using the more easily measurable T. The darkest (black) regions are 
obtained assuming that we have 10% (5%) uncertainty in R (T), the dark (green) regions 
are for an uncertainty of 20% (10%) in R (T), while the lightest (cyan) regions are for 
uncertainty of 10% (5%) in R (T). We show these regions in Fig. Elfor the four benchmark 
sets of mixing parameter cases that we have considered in this paper. The upper left hand 
panel is for TBM mixing, the upper right hand panel is for the "EXl" oscillation scenario, 
the lower left hand panel is for the "EX2" oscillation case and the lower right hand for 
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Marginalized over oscillation parameters| 





Figure 6: For each ritrue on the x-axis, the shaded region shows the values of Ufn which 
can be excluded at 2a. We assumed 5% (black/darkest), 10% (green/dark) and 20% 
(cyan/lightest) uncertainty on AT/T. The two panels are for the benchmark values TBM 
and EX2 of the oscillation parameters assumed in the data. Scenarios EXl and SPL have 
no sensitivity at all. 



the special case ("SPL") of trimaximal mixing in the second row. For T we show the 
cases for "TBM" and "EX2" only in Fig. El The way we have generated these figures is 
the following: for each value of ntme we generate the data for the particular oscillation 
parameter set. For instance for the upper left hand panel of Fig. [5] the data is always 
generated for the TBM mixing parameters. This data is then fitted with any value of 
Ufit, while all oscillation parameters are allowed to vary within their 3a allowed range. 
From Fig. O we find that the neutrino telescopes have reasonable sensitivity to the initial 
flux composition if R can be measured with reasonable precision. The sensitivity depends 
sharply on the true value of the oscillation parameters as well as on the true value of n. 
We also note that irrespective of the true oscillation parameter scenario, there is a narrow 
band around ntme ~ 1, where we have no sensitivity to the initial flux composition. From 
Fig. [6] we see that for TBM mixing and EX2, T returns a sensitivity which is not bad. In 
particular, for TBM mixing, which is the currently favored scenario, one finds reasonable 
sensitivity to the flux composition if ntme ~ 1- On the other hand we have numerically 
checked that if the true neutrino mixing turns out to be compatible with the EXl and 
SPL case, then it would be almost impossible to determine the UHE neutrino initial fla- 
vor composition using T alone. In other words, there are no shaded regions for SPL and 
EXl in the Utme—i^fu space. This is why we do not show the panels for these cases in Fig.O 

Of course it is expected that the current uncertainties on the oscillation parameters will 
decrease [31], as better neutrino oscillation experiments are built and more data is accu- 
mulated. This would improve further the prospects of measuring n at neutrino telescopes. 
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6 Conclusions and Summary 



Observation of neutrinos with very high energies coming from astrophysical sources has 
been long overdue. The observed UHE fluxes have a flavor composition on Earth which 
depends on both their initial composition at source as well as on neutrino flavor oscillations. 
Therefore, it should be possible to study both kinds of physics with these observations - 
neutrino physics as well as physics concerning the UHE neutrino sources. Unfortunately, 
since there are uncertainties in both neutrino mixing parameters as well as knowledge 
on the type of sources of these UHE neutrinos, attempts on determining one are always 
plagued by uncertainties on the other. 

In this paper we attempted to probe the potential of the future neutrino telescopes in 
deciphering the flavor composition of the UHE neutrino flux at source. We did this by 
proposing a very simple and minimal parametrization of the flavor composition at source. 
We parameterized the initial flux composition of high energy neutrinos simply as 

(*°:<:*°) = (1:^:0). (31) 

The parameter n can take any value from to oo. Specific values of n describe the pure 
known sources. The UHE neutrino flux depends on a number of astrophysical factors and 
its exact flavor composition will depend on the process which generates this flux. Neutrinos 
coming from pion decay will have a flavor ratio of (1 : 2 : 0). Likewise a muon-damped 
source will give a ratio of (0 : 1 : 0), a neutron beam source will give (1:0:0) and a 
charm source (1:1: 0). One expects impurities in these relations, i.e., they will not be in 
the exact form given above. The astrophysical source producing the UHE neutrinos could 
have a combination of some or all such processes producing these neutrinos. Also, unless 
there are close-by astrophysical sources, what will be detected at the neutrino telescope is 
most likely going to be a diffuse flux coming from a combination of different astrophysical 
sources, and hence it is expected that this observed flux would come as a combination 
of different neutrino generating decay processes. Since the number of u-,- produced at the 
source is negligible, our one parameter description is the most economical and general way 
of probing the initial flavor composition of the UHE neutrinos. 

The observed flavor ratios depend on the neutrino mixing parameters in addition to the 
property of the source. We studied two of the most used observed flux ratios T — 
^ix/i^e + + ^r) and R — $e/^r- We wrote down approximate analj^ic forms for 
these flux ratios in terms of n and the mixing parameters. Using exact numerical results 
we showed the uncertainty in T and R coming from the current uncertainty on the mixing 
parameters. It was shown that the uncertainties due to mixing parameters in minimum 
around n ~ 1 — 2. We also showed that for both T and R the uncertainty due to mixing 
for small n was more than that for large n. We also pointed out a special case where 
T = 1/3 for all values of n. This scenario corresponds to the case where the second row of 
the neutrino mixing matrix has equal entries in all its three elements and illustrates that 
the ratio T alone may not suffice to probe n. 
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We next defined a simple function and expounded the potential of the neutrino tele- 
scopes in determining n and hence the UHE neutrino fiux composition at source. Since the 
forthcoming neutrino telescopes are yet to collect any data on ultra high energy neutrino 
fiuxes, and since it is not yet known what kind of total uncertainty we would have on 
the observed fiux ratios, we performed the analysis for four benchmark points in the 
mixing parameter space and by assuming different values of errors on T and R. Those 
errors cover assumptions ranging from plausible to optimistic, and allow to compare the 
prospects of neutrino telescopes. In particular, our analysis shows what kind of statistics 
and systematics are required from neutrino telescopes in order to achieve UHE source flavor 
sensitivity. We flrst studied the case by assuming a pure pionic source. It was seen that 
in this case, once the uncertainties due to oscillation parameters were taken into account, 
T as measured in neutrino telescopes could give absolutely no information about the flux 
composition at source. The measured R could still be used to exclude certain ranges of n 
and hence certain types of sources. Finally, we performed a full scan of the n space, where 
data was generated at every value of Utrue and the potential of the data to pick the right 
source n was studied. We presented 2o" contours in the rttme ~^fit parameter space, which 
give the 2a initial flux composition sensitivity of the neutrino telescope, as a function of 

^true- 

In conclusion, neutrino telescopes will provide information on the flavor composition of the 
UHE neutrino fluxes on Earth. This can be used to study the initial flavor composition of 
these fluxes at their source. Reasonable sensitivity to deciphering the correct source flavor 
composition is expected from these experiments, despite the current uncertainties on the 
mixing parameters. With projected improvements in our understanding of the neutrino 
mixing angles, the source flavor sensitivity of the neutrino telescopes will improve. 
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